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DIRECTIONAL  WAVE  SPECTRA  USING  COSINE-SQUARED 
AND  COSINE  2S  SPREADING  FUNCTIONS 


PURPOSE  :  This  Technical  Note  presents  the  basic  concept  of  ocean  vave  direc¬ 
tional  spreading,  why  it  is  important  in  coastal  engineering,  and  how  it  is 
characterized  through  parametric  relationships.  The  COSINE- SQUARED  and  the 
COSINE  2S  spreading  functions  are  included. 

INTRODUCTION  :  The  wave  energy  at  a  point  has  an  angular  distribution  as  well 
as  a  distribution  over  a  range  of  frequencies.  This  angular  distribution  of 
wave  energy  is  termed  "directional  spreading."  Spectral  representations  which 
include  both  the  frequency  distribution  and  the  angular  spreading  of  wave 
energy  are  known  as  "directional  spectra." 

Knowledge  of  directional  spectra  regarding  their  evolution  in  shallow 
water  is  important  to  coastal  engineers.  More  accurate  vave  predictions  can 
usually  be  achieved  when  directional  spectra  are  considered.  Indications  are 
that  spectral  models  without  spreading  can  overpredict  significant  rave 
heights  during  refraction  by  as  much  as  20  percent  when  compared  to 
directional  spectral  models  (Collins,  et  al .  1981).  This  can  have  a 
significant  effect  on  coastal  structure  design  for  nondepth-limited  waves. 
Forristall,  et  al .  (1978)  concluded  that  a  proper  description  of  storm  vave 
flow  velocities  is  dependent  on  the  correct  representation  of  the  directional 
spread  of  wave  energy.  Neglecting  the  spreading  results  in  overprediction  of 
flow  velocities.  The  concept  of  directional  spreading  accounts  for  some 
important  aspects  of  wave  behavior  which  are  known  to  occur.  An  example  is 
given  in  Appendix  C  of  ETL  1110-2-305  (Vincent  and  Lockhart  1983),  where  the 
analysis  is  based,  in  part,  on  directional  spreading.  Finally,  more  realistic 
models  of  coastal  processes  are  being  developed  which  require  directional 
spectra  as  input. 
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The  ability  to  represent  directional  spectra  in  parametric  form  permits 
simulation  of  real  seas  in  numerical  rave  models,  parametric  rave  models,  and 
in  CERC's  new  directional  spectral  rave  basin. 

NATURAL  PROCESSES :  In  the  initial  stage  of  rave  generation  the  rave  crests 
are  very  short,  and  the  raves  can  be  seen  propagating  in  many  directions  dif¬ 
ferent  from  the  wind  direction.  Waves  on  a  small  lake  illustrate  this  initial 
growth  stage.  Over  longer  fetch  distances  the  raves  evolve  into  longer- 
crested  waves.  The  wind  generated  raves  develop  so  as  to  propagate  generally 
in  the  direction  of  the  wind;  but  there  is  still  an  angular  spread  of  energy 
about  the  mean  direction,  and  the  representation  of  the  rave  field  with  direc¬ 
tional  spectra  is  useful.  Waves  generated  by  a  local  storm  illustrate  this 
stage  of  rave  growth.  Finally,  the  raves  leave  the  generation  area  and  become 
swell  raves.  Ocean  swell  raves  which  have  propagated  long  distances  are  long- 
crested,  and  the  angular  spreading  of  energy  is  very  slight.  However,  repre¬ 
sentation  of  swell  in  terms  of  directional  spectra  is  still  a  useful  concept. 
Long-period  swell  arriving  on  the  West  Coast  is  a  good  example  of  long-crested 
raves. 

When  raves  approach  a  coastal  area,  other  factors  begin  to  modify  the 
directional  nature  of  the  ravefield.  For  example,  refraction  due  to  varia¬ 
tions  in  depth  can  act  to  concentrate  or  disperse  rave  energy.  Currents  can 
also  influence  the  directional  spread  of  rave  energy  in  shallow  rater  as  well 
as  in  deep  rater. 

PARAMETRIC  SPREADING  FUNCTIONS:  The  basic  premise  of  all  parametric  spreading 
functions  is  that  the  single-peaked  directional  spectrum  is  adequately 
described  by  the  product  of  two  functions: 

E(f,  0)  =  E(  f  )*  D(  f  ,  0) 

where 

E(f,  0)  »  directional  spectral  density  function 

E(f)  «  one-dimensional  energy  spectral  density  function 
D(f,  0)  =  angular  spreading  function 

f  «  frequency  in  hertz 
0  =  direction  in  radians 

The  formulation  of  D(f,  0)  requires  that 
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E(f)  D(f,  0)  d  0df  =  J  E(f)  df 

0 

so  that  the  total  energy  in  the  directional  spectrum  must  be  the  same  as  the 
total  energy  in  the  corresponding  one-dimensional  spectrum.  This  parameter¬ 
ization  can  effectively  represent  the  directional  nature  of  a  wavefield  in  the 
absence  of  complicating  influences  such  as  a  large  change  in  wind  direction  or 
the  propagation  of  swell  into  the  generation  area,  which  results  in  a  bimodal 
spectrum.  Several  idealized  spreading  functions  are  outlined  below. 


1 .  Co  s  ine-  Squa  re  d 

cos2(e-en)  ,  for  c-  ~  +  en)  <  0  <  (■£  +  0.) 

D(0)  4’  0  2  0  2  0 

0  ,  otherwise 

where  0  =  mean  wave  direction  in  radians, 

o 

The  cosine-squared  formulation  is  extremely  simple  because  it  is  neither  a 
function  of  frequency  nor  widespeed.  It  can  be  used  to  parameterize  the 
directional  spreading  of  wind  seas.  Similar  formulations  can  be  derived 
(cosine-fourth,  for  example)  by  changing  the  value  of  the  exponent  and 
adjusting  the  coefficient  so  that 


f 


D(0)  d  0  =  1 


Most  CERC  spectral  wave  models,  including  the  models  used  by  the  Wave  Infor¬ 
mation  Study  (WIS)  (CETN-I-19),  use  the  cosine-fourth  spreading  function  to 
represent  locally  generated  waves  (seas). 

2.  Cosine-2S 


D(f,  0) 


(2s  - 


TT 


r2(s  +  l) 
r  (2s  +  l) 


3 


where 

T  =  the  Gamma  function 

0  =  the  mean  wind  direction 

o 

s  =  the  spreading  parameter  which  is  a  function  of  frequency  and 
wind speed 

This  spreading  function  was  proposed  by  Longue t-Higgins ,  et  al  •  (1963).  Three 
of  the  most  commonly  used  and  referenced  formulations  for  the  parameter  s 
are  as  follows:  (a)  s  =  constant;  (b)  Mitsuyasu,  et  al .  (1975);  and  (c) 
Hasselmann,  et  al.  (1980). 

The  first  parameter,  (a),  is  comparable  to  consine-nth,  which  is  a 
generalization  of  the  cosine-squared  spreading  function  presented  earlier. 

The  change  in  directional  spread  as  a  function  of  s  is  illustrated  in 
Figure  1.  Increasing  the  power  of  the  cosine  causes  a  narrowing  of  the 
directional  spread.  Swell  is  more  typically  represented  by  the  narrow  spreads 


v(e) 


Figure  1.  Idealized  angular  distribution 
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D( 0)  =  G(s)  cos 
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in  Figure  1.  Most  numerical  and  physical  modeling  of  directional  spectra  make  use  of  this 
formulation.  A  simple  approximation  to  this  formulation  is  discussed  in  CETN-I-6(1985). 

The  last  two  formulations,  (b)  and  (c),  are  based  on  field  data  in  deep  water.  Both  include  windspeed  and 
peak  spectral  frequency  in  the  parameterization,  although  windspeed  is  not  a  critical  part  of 
Hasselmann's,  and  it  could  be  left  out.  While  these  latter  two  parameterized  spreading  functions  more 
accurately  depict  nature,  they  generally  are  not  used  in  numerical  models  because  they  are  inefficient 
numerically;  and  their  increased  accuracy  is  lost  in  the  coarse  angular  resolution  of  the  models. 

APPLICATION:  Practical  applications  of  parametric  directional  spreading  functions  include  numerical 
and  physical  modeling  of  waves  and  analysis  of  directional  field  data.  Field  measurement  analysis  often 
involves  selecting  a  value  for  n  or  s  in  the  parameterized  spreading  function  which  optimizes  the  fit  to 
measured  directional  attributes.  When  measurements  of  directional  spectra  are  unavailable  or 
impractical,  reasonable  approximations  for  numerical  or  physical  modeling  can  be  made  by  combining  a 
spreading  function  with  either  measured  one-dimensional  spectra  or  with  empirical  one-dimensional 
spectral  formulations,  such  as  the  JONSWAP  spectrum  in  deep  water  or  the  TMA  spectrum  in  shallow 
water  (Hughes  1984).  The  composite  directional  spectra  then  provide  useful  input  for  shallow- water 
spectral  transformation  numerical  or  physical  models.  Present  shallow-water  numerical  models  allow 
calculation  of  directional  spectra  in  the  nearshore  region;  and  they  account  for  effects  such  as  refraction, 
bottom  friction,  nonlinear  interactions,  wave  breaking,  and  currents. 

ADDITIONAL  INFORMATION.  For  additional  information  contact  Ms.  Jackie 
Pettway,  Jackie.W.Pettwav@usace.army.mil  of  the  Coastal  and  Hydraulics  Laboratory. 
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